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Abstract 

We derive raising and lowering operators for orthogonal polynomials on the unit circle 
and find second order differential and q-difference equations for these polynomials. A general 
functional equation is found which allows one to relate the zeros of the orthogonal polynomials 
to the stationary values of an explicit quasi-energy and implies recurrences on the orthogonal 
polynomial coefficients. We also evaluate the discriminants and quantized discriminants of 
polynomials orthogonal on the unit circle. 
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1. Introduction. Let w(z) be a weight function supported on a subset of the unit circle and 
assume that w is normalized by 



(1-1) / w(C)^ = l. 



Let <t> n {z) be the polynomials orthonormal with respect to w(z), that is 

d( 



(1.2) / <MC)0n(CMC)v = < W- 
J\(\=i < 

A general background to orthogonal polynomial systems defined on the unit circle can be found 
in the monographs [33], [|l(|, and ||, while more recent surveys are to be found in [11], [12] and 
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p4| ] and from an interesting perspective, in the course notes of [23|. In this work we first derive 
raising and lowering operators for <j) n under certain smoothness conditions on the weight function 
then use these first order operators to derive a linear second order differential equation satisfied 
by the orthogonal polynomials. This will be done in §2. These results are unit circle analogues of 
the results of Bauldry Q , Bonan and Clark Q , and Chen and Ismail || . The external field 119] , 
[p7|1 is the function v defined by 



(1.3) w(z) = exp(— v(z)). 



We illustrate these general results by three examples - the circular Jacobi polynomials, the Szego 
polynomials and the orthogonal polynomial system defined by the modified Bessel function. 

Flowing from the results in §2 we derive a functional equation and relate this to the zeros of 
orthogonal polynomials defined on the unit circle in §3. This is the analogue of the electrostatic 
interpretation of the zeros of orthogonal polynomials defined on the real line, but in this case 
analyticity means that the quasi-energy function derived has stationary points at the zeros which 
are saddle-points, not minima. This functional equation implies a general relationship on the or- 
thogonal polynomial system, which is usually expressed as a recurrence relation on the polynomial 
coefficients. 

In §4 we derive g-analogues of §2. The external field is now the function u defined through 



(1.4) (D q w)(z) = -u(qz)w(qz), 
where D q is the g-difference operator 

(1.5) (D, m := ' iZ \-J£ Z) . 
In other words 



(1.6) w{z) = w(qz)[l — (1 — q)zu(qz)], \z\ = 1. 



Recall that the discriminant D(f n ) of a polynomials f n is defined by 



(1.7) D{f n )= 1 2n - 2 J] ( Zj -z k ) 2 , if f n (z)= 1 H(z-z j ). 

l<j<k<n j=l 

Stieltjes [30], [31] and Hilbert [15] evaluated the discriminants of the classical orthogonal poly- 
nomials of Hermite, Laguerre, and Jacobi. Schur |28|] gave an interesting lemma about general 
orthogonal polynomials which implies the Stieltjes-Hilbert results when applied to the Hermite, 
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Laguerre, and Jacobi polynomials. In §4 we prove an analogue of Schur's lemma for polynomials 
orthogonal on the unit circle and use it to give a general theorem on the evaluation of discrimi- 
nants of orthogonal polynomials on the unit circle. This is the unit circle analogue of the results 



in Jlq ]. It was observed in [21] that, in general, the discriminant (1.7) of g-orthogonal polynomials 
does not have a closed form. The appropriate discriminant for discrete q-orthogonal polynomials 
is 

(1.8) D(f n , q) = 7 2 - 2 q^' 2 H (q^ Zj - z^ 2 )^ 2 Zj - z k q^ 2 ), 

l<j<k<n 

if /„ is as in (|L7]). The above discriminant also has the alternate representation 

(1.9) D(f n , g) = 7 2 ™- V (n " 1)/2 II + + 

l<j<k<n 

In particular for a quadratic polynomial Az 2 + Bz + C the g-discriminant is qB 2 — (1 + q) 2 AC. 

In §5 we give an expression for the ^-discriminant of polynomials orthogonal on the unit 
circle in terms of the coefficients in the recurrence relations satisfied by the polynomials. As an 
illustration we evaluate the g-discriminant of the Rogers-Szego polynomials [|32||, 



2. Differential Equations. Recall that if / is a polynomial of degree n then the reciprocal 
polynomial is 

n n 

(2.1) r(z):=J2^z n ~ k , if f{z) = Y j a k z\ and a n + 0, 

k=0 k=0 



pq|. Let (p n satisfy (1.2) and 

(2.2) <t>n{z) = K n z n + l n z n ~ x + lower order terms, K n > for n > 0. 
Then the n 's satisfy the recurrence relations f33] , (11.4.6), (11.4.7)] 

(2.3) K n Z(j) n (z) = K n+ l4> n+ l{z) - (j) n+1 (0)(j)* n+1 ( Z ), 

(2.4) K n (j) n+1 (z) = K n+ lZ(j) n ( Z ) + (j) n+ i(0)(j)* n ( Z ). 



If we eliminate </>* between ( |2.3] ) and (2.4) we get the three term recurrence relation (XI. 4, p. 91 
in 111 



(2.5) K n (j) n (0)(j) n+1 (z) + K n -l(f) n+ l(0)z(j) n -l(z) = [K n (p n+1 (0) + (0)z] n (z) 
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The re's and <f) n (0) are related through |33], (11.3.6)] 

n 

(2.6) k\ = £ | 



^0)| 2 . 



fc=0 



Thus re n (> 0) can be found from the knowledge of |0fc(O)|. By equating coefficients of z n in (|2.5D 
and in view of (12.21) we find 



(0) + re n _ 1 n+ i(O) = re n n+ i(O) + (0). 

Thus 



(2.7) re„/ n+ i = K n+ il n + n (O)^ n+ i(O). 

From (|2.7|) it is possible to express in terms of the re's and </>j(0)'s 
1 PM+i (°) 



(2.8) Z n = re n 



The analogue of the Christoffel-Darboux formula is 



(2.9 > fawfaiz = : = ■ 

^— ' 1 — az 

k=0 



Theorem 2.1 Let w(z) be differentiable in a neighborhood of the unit circle, has moments of all 
integral orders and assume that the integrals 



ld=i Z ~C < 

exist for all integers n. Then the corresponding orthonormal polynomials satisfy the differential 
relation 

(2.10) <j>' n {z) = n^^_i(*)-<(*) / V ' {z) ~ V ' i0 M0'MC~MC)d( 



v'(z)-v'(0 



+i(t>n{z) / W ^> n (C)^n(CMCK 

'ic[=i «-C 
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Proof. Using the orthogonality relation (|1.2|) we express 4>' n {z) in terms of the </>fc's as 



n-l 



H h\=i < 



k=0 



t/(C)0fc(C) + C<MC) + CV fc (C)l MCMC)^, 

J it, 



where we have integrated by parts, then rewritten the derivative of the conjugated polynomial in 
the following way 



d 



-2- 



(2.11) -0 n ( C ) = -c <(C), 



and used the fact C = VC f° r I CI = 1- Now the orthogonality relation fll.2j ) and (2^2) give 



4>M 



n-l 



Kl=i < 

+^n-l(z) f [C^-l(C) +C 2 <-l(C)l 4>n(CKC) 

7|ci=i L 
t/(*)-t/(0 



^n— 1 



0n(C) faiOfti*) - MOM*) w (0 



+<j) n -i(z) 



+ (n-l)- 



"n— 1 



This establishes ( p.lOj ) and completes the proof. 

We next apply (pit ) to eliminate </>* from ( 2.1C|) , cissuniing 0n(O) 7^ ^* The result is 

v'(z)-v'(C) 



(2.12) 



Kn 4>n{0) J\C\=1 Z-Q 



-i<t>n{z) 



v'(z)-v'(C) 
Kl=i 



</>n(C) 



tfn(C) " T^n(C) 



Observe that <p n {C) ~ ^n(C) is a polynomial of degree n — 1. 
Let 

- t/(C) 



(2.13) A n {z) 

(2.14) S n (z) 



f^n—l . - ^n— 1 
«™ <^n(0) 7|C|=1 z - C 



t/(f) ~ v'tt) 

ld=i «-C 



0n(C) 



0„(C)#i(CMCK, 



4>n{0) 
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For future reference we note that Aq = Bq = and 

(2.15) A 1 {z) = n 1 -<j )1 {z)v'{z)-^-M x (z), 

(2.16) B 1 (z) = - v '( z )-^ Ml (zl 

<Pi (0) 

where the first moment Mi is defined by 

7ld=i < 



Now rewrite ( j2.12 ) in the form 



(2.18) (f>' n (z) = AnWfa-^z) - B n {z)<j) n {z). 



Define differential operators L n ^\ and L n ^ by 
(2.19) L nil = ^- + B n (z), 

, n on ^ r d o / \ , An-l^^n-l . A n -i(z)K n l/> w -i(0) 

(2.2U) L n2 = -j B n -i(z)-\ 1 — — . 

dz ZK n _ 2 K n _ 2 0n(O) 



After the elimination of 4> n -\ between ( [2.12 ) and (|2.5[) we find that the operators L n ^\ and L H)2 



are annihilation and creation operators in the sense that they satisfy 

(2.21) L n ^ n (z) = Mz^iiz), L n ^ n ^{z) = ^M. ^f )K "-V „(4 
Hence we have established the second order differential equation 

(2.22) L n2 — — -L„,i n (z) = — (p n {z), 

\A n (z) ) Z 4> n (0)K n -2 

which will also be written in the following way 

(2.23) ^ + P{z)^ n + Q{z)<t> n = Q. 

It is worth mentioning that, unlike for polynomials orthogonal on the line, L* 1 is not related 
to L n 2- In fact if we let 



(2.24) (/,<?):=/ /(C)ff(CMC)^, 
J\C\=i s 
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then in the Hilbert space endowed with this inner product, the adjoint of L n i is 



(2.25) (L* nl f){z) = z 2 f(z) + zf(z) + [v(z) + B n (z)]f{z). 



To see this use integration by parts and the fact that for |£| = 1, g(() = g(l/(). 

Example 1 The circular Jacobi orthogonal polynomials (CJ) are defined with respect to the 
weight function 

for real a appropriately restricted. We find these to be classical in the sense of being related to 
classical orthogonal polynomials defined on the real line and therefore possessing their properties. 
They arise in a class of random unitary matrix ensembles, the CUE, where the parameter a is 
related to the charge of an impurity fixed at z = 1 in a system of unit charges located on the unit 
circle at the complex values given by the eigenvalues of a member of this matrix ensemble |3f| • 
The orthonormal polynomials are 

(2.27) <t> n (z) = - r =£k== 2 F 1 (-n,a+l-,-n+l-a-z) , 

v / n!(2a+l) n 

and the coefficients are 

(2.28) Kn = ; (Q+1)n n>0, 

Vn!(2a + l) n 

na 

(2.29) l n = K n n > 1, 

n + a 

(2.30) 4> n (0) = -^—K n n>0. 

n + a 

The reciprocal polynomials are 

(2.31) </>*(z) = -= jy=^L= 2 F 1 (-n,a;-n-a;z). 

^n\(2a+l) n 

Using the differentiation formula and some contiguous relations for the hypergeometric functions, 
combined in the form 

. d . , n(n + 2a) „ , 

2.32 {l-z)— 2 F 1 {-n,a+l;l-n-a;z) = '- 2 F X 1-n, o+l; 2-n-o; z) 

az n — l+a 

— n 2 F\{—n, a + l;l — n — a; z), 

one finds the differential-recurrence relation 
(2.33) (1 - z)(j)' n = -n^ n + [n(n+2a)] 1 / 2 0n-i, 
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and the coefficient functions 



. / \ \/n(n+2a) 

(2.34) A n (z) = V ±— \ 

ft 

(2.35) B n {z) = ^— . 

The second order differential equation becomes 

^oqk\ / // I / / fl-n-a 2a+l\ n(a+l) 

(2.36) *„ + *„| — = 0. 

Example 2. We consider a generalization of the previous example, to the situation where 

(2.37) w(z) = 2 - l - 2a - 2b , r( ^^ 1} ; J l - z\ 2a \l + z| 2fe , 

v ; w r(o+i/2)r(6+i/2) 1 11 1 

with x = cos 6, and the associated orthogonal polynomials are known as Szego polynomials ]3 
They are related to the Jacobi polynomials via the projective mapping of the unit circle onto the 
interval [—1, 1], z t— > -[z+z^ 1 ] = x = cos 6, 

(2.38) z"> 2n (z) = AP(-V^-V2)(1 [z+z _ 1]} + l B[z _ z _ 1]p (^ 1Ab+1/2)( l [z+z _ 1]}j 

(2.39) ^^(z) = CP^-V^-^) + Id[.-,-1]p(^ A6+1/2) (^-+^ 1 ])- 

In their study of the equilibrium positions of charges confined to the unit circle subject to log- 
arithmic repulsion Forrester and Rogers considered orthogonal polynomials defined on x which 



are just the first term of 2.38 . Using the normalization amongst the even and odd sequences of 
polynomials, orthogonality between these two sequences and the requirement that the coefficient 
of z~ n on the right-hand side of ( p. 39 ) must vanish, one finds explicitly that the coefficients are 



(2.40) A 



n\(a+b+l) 



(o+l/2) n (6+l/2) n 



1/2 



(2.41) B = -A, 



2 



(2.42) C = n 



(2.43) D = n+a+b C . 



(n-l)!(q+6+l) n _i 
(a+l/2) n (6+l/2) n 

2n 



1/2 
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Furthermore the following coefficients of the polynomials are found to be 

o-2n (a+b+l) 2n 

v / n!(a+5+l)„(a+l/2) n (6+l/2) n ' 

2 l-2n (a + 6+l) 2w -l 

v / (n-l)!(a+6+l) n _i(a+l/2)„(6+l/2)„' 
2n 2^6 K2 "' 

(2re " 1) 2n+7+6-l K2 "- 1 ' 
a + b 



(2.44) 


K2n 


(2.45) 


^2n— 1 


(2.46) 


^2n 


(2.47) 


^2n-l 


(2.48) 


02n(O) 


(2.49) 


02n-l(O) 



2n+a+b 
a — b 

«2n-l- 



2n+a+6-l 
The three term recurrences are then 



(2.50) 2(a-6)v / n(n+a+6)02n(z) + 2(a+6) y/(n+a-l/2)(n+b-l/2)z<f>2 n - 2 (z) 

= [{a+b){2n+a+b-l) + (a-6)(2n + a+6)z]0 2 n-i(^), 

and 

(2.51) 2(a+b) y / (n+a-l/2)(n+b-l/2)(f> 2n _ l (z) + 2(a-b)V("-l)(ra+a+6-l)z02n-3(2) 

= [(a-b)(2n + a+b-2) + (a + 6)(2n + a+6- l)z]<fen_ 2 (2), 

when a ^ b and both these degenerate to (f>2n-i{z) = z$2n-i{z) when a = b. For reference the 
reciprocal polynomials are 

(2.52) z-^Uz) = AP(-V^-V2) (![*+*"!]) - ^[z-z-^P^'^cilz+z- 1 ]), 

(2.53) s^-iM = ^^/^-^(l _ ^_ z -i]^iAW/2) ( l [z+z -i ]} . 

Using the differential and recurrence relations for the Jacobi polynomials directly one can find 
the appropriate coefficient functions for the Szego polynomials to be 



(2.54) 


A 2n -i{z) = 


(2.55) 


B 2n -i(z) = 


(2.56) 


A 2n {z) = 


(2.57) 


B 2 n(z) = 



„ r, i—n ; r--a — b + (a + b)z 

2V(n+a-l/2)(n+6-l/2) (a _ fc)( \_^ ) , 

Aab + (2n-l) [a+b + (a-b)z\ 



= 2n 



(a-b)(l-z 2 ) 

7T a + 

-a+b — 
(a- 

a— 6 + (a+6)z 
(a+fe)(l-^ 2 ) ' 



^— ; rra+b + (a—b)z 

2Vn(n + a + b) (a + ft)(1 _ z2) , 
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again when a 7^ 6. 

Example 3. Consider the weight function 



(2.58) 



w[z 



^exp^z + z- 1 ]), 



where I v is a modified Bessel function. This system of orthogonal polynomials has arisen from 
studies of the length of longest increasing subsequences of random words || and matrix models 
[|26}|, [16]. In terms of the leading coefficient one has the Toeplitz determinant form 

(2-59) 4W = Io(t)^^^. 

aet{lj-k{t))o<j,k<n 

The first few members of this sequence are 



(2.60) 
(2.61) 
(2.62) 
(2.63) 



l(0) 



Kl 



hit) 

io(ty 



«2 



i(0) 



I (t)(p (t)-If(t)) 



(/„(*) - I 2 {t)) [/ (t)(/ (0 + h(t)) - 2lf (*)] 
I (t)I 2 (t)-I*(t) 



m-i 2 (t) ■ 

Gessel [jDJ has found the exact power series expansions in t for the first three determinants which 
appear in the above coefficients. Some recurrence relations for the corresponding coefficients of 



K2 



the monic version of these orthogonal polynomials have been known |26j, [16], |3J] and we derive 
the equivalent results for K n , etc. 



Lemma 2.2 ( pq] ) The reflection coefficient r n (t) = ^ n (0)/K n for the modified Bessel orthogonal 
polynomial system satisfies a form of the discrete Painleve II equation, namely the recurrence 
relation 



(2.64) 



n r r , 
2— 



t 1 



r n +i + r n -\, 



forn>\ and r (t) = I, n(t) = -h(t)/I (t). 



Proof. Firstly we make a slight redefinition of the external field w(z) = exp(— v (z + 1/z)) for 
convenience. Employing integration by parts we evaluate 



- J V'(C + 1/0(1 -l/C 2 )^n+l (C)tfn(CMC)- 



0n+l(C)CVn(C) + </>n+l(C)C<M0 - 0n+l(C)0n(C)l ™(C) ^ 



(2.65) 



(n+1) 
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for general external fields v(z) using Ql.2| ) and ( |2.2| ) in a similar way to the proof of Theorem 2.1. 
However in this case v'(( + 1/Q = —t/2, a direct evaluation of the left-hand side yields 



(2.66) 



1 

-t 



In 



n lj n+2 



Kn+1 Kn+lKn+2 



and simplification of this equality in terms of the defined ratio and use of ( |2.8| ) gives the above 
result. 

There is also a differential relation satisfied by these coefficient functions or equivalently a 
differential relation in t for the orthogonal polynomials themselves jl6| , pi] . 



Lemma 2.3 The modified Bessel orthogonal polynomials satisfy the differential relation 



(2.67) 2-6 n {z) 
dt 



h(t) <f> n+1 (0) 



I (t) K n+1 <p n (0) 



<Pn(z) 



K n—1 



-12 



for n > 1 and ifoiz) = 0. The differential equations for the coefficients are 



(2.68) 
(2.69) 



2 d 
2 d 



<Mo) dt 



(0) 



h(t) n+ i(O)0 n (O) 
Ji(t) , 0„+i(O) (0) 



+ 



(0) 0n(O) 



for n > 1 . 



Proof. Differentiating the orthonormality relation (|l.2j ) with respect to i one finds from the 
orthogonality principle for m < n — 2 that 

(2.70) + ^ Z( Pn( z ) = a n <j) n+ i(z) + b n (p n (z) + c n ^ n _i(z) 

for some coefficients a n ,b n ,c n . The first coefficient is immediately found to be a n = |«„/« n+ i. 
Consideration of the differentiated orthonormality relation for m = n — 1 sets another coefficient, 
c n = — ^K n _i/K n , while the case of m = n leads to 6 n = ili(f)/Io(i). Finally use of the three- 
term recurrence ( glj ) allows one to eliminate (/> n+ i(z) in favor of 4> n {z), 4> n -i{ z ) an d one arrives at 
Q2.67|) . The differential equations for the coefficients K n ,(j> n (0) in ( p.68 , p.69|) follow from reading 
off the appropriate terms of fl2.67| ). 

Use of the recurrence relation and the differential relations will allow us to find a differential 
equation for the coefficients, and thus another characterization of the coefficients. 
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Lemma 2.4 The reflection coefficient r n {t) satisfies the following second order differential equa- 
tion 



— 

^ ' ' dfi r ' 1 2 V r n + 1 ' r n - 1 J V dt 



1/1 1 

+ 



\ f d \ Id 2 n 2 r n 

){dt rn ) --tdt rn - rn{1 - r - ) + w— 



with the boundary conditions determined by the expansion 
for n > 1. T/ie coefficient r n is related by 

to z n (t) which satisfies the Painleve transcendent P-V equation with the parameters 



n 



(2.74) a = -/3 = —, 7 = 0, 5 = -2. 



Proof. Subtracting the relations ( 2.6§| ,2.69) leads to the simplified expression 



2 d 

(2.75) r n+ i - r n _i = K — ry, 



\-r\dt 



which should be compared to the recurrence relation, in a similar form 



(2.76) r n+1 + r n _! = -2 



n r r , 



tl-rl 



The differential equation ( p.71[) is found by combining these latter two equations and the identi- 
fication with the P-V can be easily verified. 

As a consequence of the above we find that the coefficients for the modified Bessel orthogonal 



polynomials can be determined by the Toeplitz determinant ( p.59|) , by the recurrence relations 
( [2.76 ) or by the differential equation ( |2.7l| ). An example of the use of this last method we note 



(2.77) = I (i) [1 - r 2 n (t)]- 1/2 e W f-n J ~ 



ds rl(s) 



We now indicate how to find the coefficients of the differential relations, A n (z), B n {z) and 
observe that 



v'(z) - v'(Q 
z-C 



1 1 

<2 
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The above relationship and ( p. 12 ) yield 



(2.78) <f>' n (z) 



071-1 C?) 



n + 



2z 



2 0„(O) 



IC]=i 



2z 



/ MOC 



Kl=i 



MO 



K(0 



w(C) 



K' 



Easy calculations using Q2.2|) give 



and 



MC) - ==M0 



MO = 0n+l(C) , 



K n -1 0n-l(O) 



»(0) 



«n n (O) 
^rJ"n ln+1 



4> n (0 + lower order terms. 



kn+l K n 



4>n{0 + lower order terms. 



These identities together with ( 2.7S| ) establish the differential-difference relation 
(2.79) <p' n (z) 



Kn—l 



t tK n _i0 n _i(O) t (j) n+ l(0)4> n (0) 

n + — + -■ 



2z 2 K n (p n (0) 2 K n+i K n 



0n-l(z) 



t K w _i n -l(O) 
2z K n n (O) 



3. Functional Equation and Zeros. In this section we continue the development of the 
previous discussion of the differential relations satisfied by orthogonal polynomials to find a func- 
tional equation and its relationship to the zeros of the polynomials. Expressing the second order 
differential equation ( 2.22| ) in terms of the coefficient functions A n (z) and B n (z) we have 



(n t\ J." i J D i D A' I A K n-1 Ai-l K n </>n-l(0) . \ , 

(3.1) (j) n + <B n + - A n /A n — — A n _i ) cj) n 

+ < B' n — B n A' n /A n + B n B n ^i — 

L K n- 



K n— 1 A n —\B n 

-2 Z 

K n n _l(O) K n _i 0„-l(O) A n _ X A n 

-A n -iB n H — - — S = 0. 



(0) <£n(0) JZ 

Now by analogy with the orthogonal polynomials defined on the real line the coefficient of the <p' n 
term above can be simplified. 

Theorem 3.1 Given that v(z) is an meromorphic function in the unit disk then the following 
functional equation holds 

; n\ did K n-1 A n —\ K n n _l(O) . . . _j . 

(3.2) B n + B n _i A n _i = -(n - l)z - v (z). 

K n „ 2 2; K n -2 <Pn(°) 
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Proof. From the definitions ( |2.13 , 2Tl4| ) we start with the following expression 

(0) 



B n + B 



n-l 



K n _ 2 Z 

-(n-l) 



n-l 



1 K n n -l(O) 



+ t 



2-C 



0n(O) 



y n<^ n - <Pn-Wn-l 
. K 



0n(O) 



1 n 



0„(O) 



Employing the recurrences ( |2.4|J2.3| ) , and the relation amongst coefficients ( |2,6[ ) one can show that 
the factor in the first integral on the right-hand side above is 



<>n + 



,(0) 



J nYn 



C,4>n-l(t>* n -l = -4>n4>n + 4>„4>n- 



Now since \C\ 2 = 1, one can show that the right-hand side of the above is zero from the Christoffel- 
Darboux sum (|2.9|). Consequently our right-hand side is now 



1 K n 0n-l(O) 



-(n-l) 



z K n— 1 



»(0) 



Taking the first integral in this expression and using the recurrence (2.4) and the decomposition 
(4> n -i = K n -\/ K n 4> n +i: n -i where ir n £ H n , II n being the space of polynomials of degree at most n, 
we find it reduces to —i(ft n (0) / n n from the normality of the orthogonal polynomials. Considering 
now the second integral above we integrate by parts and are left with 



<f>'n-l<t>n-MOdC + / 0n-l<A„_XCR, 



and the first term here must vanish as 4> n _i can be expressed in terms of cf> n -i,(f) n from (2.4) 
but <p' n _i £ n„_2- The remaining integral, the second one above, can be treated in the following 
way. Firstly express the conjugate polynomial in terms of the polynomial itself via ( |2.3D and 
employ the relation for its derivative ( [2.11 ). Further noting that C0n— l = ( n ~ l)<^n— l + 7r n-2 5 
(cpn~2 = K n _ 2 /K n -i0n-i + 7T n __ 2 , and C 2 </> n -2 = ( n ~ 2)K n _ 2 /t n -i0n-i + ^n~2 along with the 
orthonormality relation, the final integral is nothing but —i(n — l)<^ n _i(0) / K n -i. Combining all 
this the final result is (|3 



Remark 1. The zeros of the orthogonal polynomial 4> n (z) are denoted by {zj}i<j< n , and are 
confined to the convex Hull of the support of the measure, namely to be strictly confined within 
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the unit circle \z\ < 1. One can construct a real function \T(z\, . . . , z n )\ from 

n —viz-) 

(3.3) T^,... >Zn ) = JJ z -»+ii-_L (z,-z fc ) 2 , 

such that the zeros are given by the stationary points of this function. One might also interpret 
this function as a total energy function for n mobile unit charges in the unit disk interacting with 
a one-body confining potential, v(z) + lnA n (z), an attractive logarithmic potential with a charge 
n— 1 at the origin, (n— 1) lnz, and repulsive logarithmic two-body potentials, — hxfe—Zj), between 
pairs of charges. However all the stationary points are saddle-points, a natural consequence of 
analyticity in the unit disk. That such this function exhibits stationary properties at the zeros can 
be seen by considering the second order differential equation which in view of the above theorem 
has the coefficient P(z), namely 

(3.4) P(z) = -{n- l)z~ l - v'(z) - A'JA n . 

This function is a perfect differential and consequently the one-body potential can be constructed 
from its integral, via the Stieltjes argument. Or alternatively one can show that the conditions 
for the stationary points of function T(z\, . . . , z n ) above lead to a system of equations 

(3.5) -„' fe )-|M_^ + 2 V ^_ = „ , = !,... ,„. 

Then the pairwise sum can be represented in terms of the polynomial f(z) = Ylj = i(z — zj) thus 

i<k<n,k& Zj ~ Zk 
and we have the n conditions expressed as 

(3-7) /'^ + {-^-^)-^||}/'N=0, Vj = l,...,n. 
The result then follows. 

Remark 2. The functional equation ( |3.2| ) actually implies a very general recurrence relation on 
the orthogonal system coefficients K n ,(j) n (0). In general if it is possible to relate the differential 
recurrence coefficients A n ,B n to these polynomial coefficients, then the functional equation dic- 
tates that equality holds for all z, and thus for independent terms in z. For rational functions 
this can be applied to the coefficients of monomials in z. 
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Remark 3. There is another way of deriving the functional equation ([03) which we now describe. 



Equation ( 2.22 ) is one way of expressing the second order differential equation for the orthogonal 



polynomials, however one can perform the elimination in the opposite order and find 

(3.8) L n+11 ( Z L n+12 ) (j) n (z) = Kn ^°\ A n+1 {z) (f) n (z), 
\A n (z) ) K n _i0 n+ i(O) 

which written out in full is 

(3.9 (j> n + i B n+1 + B n - AJA n j-rA n + -)cj) n 

1 K n _i Z K n _l0 n+ l(O) Z) 

+ \ B' n - B n A' n /A n + B n+ \B n — ■ 



Z 



Kn+1 4>n{0) . D «n 0n(O) A n A n+ i B n 

-A n B n+ l H ; — 1 ; — } (p n = 0. 



Kn-1 0n+l(O) Kn-10n+l(O) 2 2 «n-l 0n+l(O) 2 

Given that the coefficient P(z) for these two forms (|3~T,3.9) must be identical we have an inho- 



mogeneous first order difference equation, whose solution is 

(3.10) B n + £ n _! - ^Hzl^zl _ ^ ^-^ Q) A n _! = -(n - l)*" 1 + function of z only. 

This function can be simply evaluated by setting n = 1, evaluating the integrals after noting 
Bq = and the cancellations, and yields the result —v'(z). 

Example 1. We can verify that the general form for the T- function is correct in the case of the 
circular Jacobi polynomials by a direct evaluation 

n 

(3.H) t( Zi ,... ,z n ) = i\z]- n - a (i-z J r +i (z J -ir n (*i-**) 2 > 

j=l t<j<k<n 

where we have used the identity 

(3.12) |1 - z\ 2a = (1 - z) a {l - l/z) a = z- a (l - z) a {z - l) a , 

on \z\ = 1 to suitably construct a locally analytic weight function. One can show that the 
stationary points for this problem are the solution to the set of equations 



1 — n — a 2a + 1 v-^ 1 

(3.13) ,^_ + 2 ^^— = 0, l<j<n, 

Z 3 1 Z 3 Z j Z k 

so that the polynomial f{z) = nj=i( z ~~ z j) satisfies the relations 

(3.14) rfe)+/'fe){^-|^j=0. 

I z j 1 z j J 
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Consequently we find that 



(3.15) z(l - z)f"[z) + f'(z) {(l-n-o)(l - z) - (2a+l)z} + Qf(z) = 0, 



for some constant Q independent of z, but possibly dependent on n and a and is identical to the 



second order ODE ( 2.36 ). 

Example 2. Using the expressions one can verify that the identity ( |3.2| ) holds and in 

particular becomes 

/ oie s t-j . t-j K n _i^„_i K n n _i(O) ra-1 a+6 2a 26 
[6. Lb) t> n + B n -i , —\ A n-l = r 



K n _ 2 Z K n -2 4>n(0) 



1-z l+z" 



for both the odd and even sequences. Consequently the coefficients in the second order differential 
equation are 



(3.17) P n (z) 

(3.18) Q 2n (z) 

(3.19) Q 2n ~i(z) 



a + b 



2n 



n+a+b-l 2a+l 2b+l 

z 1 — z l + z a + b+(a + b)z' 

a(a + l)(l + z) 2 -b(b+l)(l-z) 2 



z(l-z 2 )[a+b + (a-b)z\ 
(2n - l)[a(q+l)(l + z) 2 + b(b+l){l-z) 2 - 2ab(l-z 2 )} + 4o6 
z(l-z 2 )[a-b + {a+b)z\ ' 



Similarly we can verify that the general form for the T-function is correct in the case of the 
Szego polynomials by using the identity 



(3.20) |1 - z\ 2a \l + z\ 2b = z- a - b (l - z) a (z - l) a (l + z) 2b , 
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to suitably analytically continue the weight function. The stationary points for this problem are 
the solution to the following set of equations 



1-n-a-b 2a + l 2b+l 
(3.21) rr^ + 



a + b 



1 — Zj 1 + Zj a ± b + (a + b)zj t rf Zj — Zk 

such that the polynomial f(z) = Y\™ = i(z — Zj) satisfies the relations 

a + b 



+ 2 V — = 0, l<j<n, 



(3.22) / ^) + ^)^l=2=£^_^±l + 2ft + l 



1 — z j 1 + Zj a db b + (a + b)zj 



0. 



Finally we find that 



..,,,■>, , , ft xjl-n-a-6 2a+l 26+1 
(3.23) / (z) + f (z) { h 



a + b 



1-z l + z a + b + (a + b)z 



+ Q(z)f(z) = 0, 
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for some constant Q independent of z, but possibly dependent on n and a. The coefficient of the 
first derivative term is identical to the expression for P(z) in (pTD . 

Example 3. One can also verify the functional relation ( |3.2| ) for the modified Bessel orthogonal 
polynomials. Forming the left-hand side of this identity we find this reduces to 

(3.24) Bn + Bn _ 1 -^^l-^L^i 4 _ 1 



(0) , t^-i(0) 
2 2z 2 ln j Kn _i n (O) 2 n (O) 2 

Now the last three terms on the right-hand side of the above equation simplify to t/2 using the 
recurrence relation (2.64), showing that the general functional relation holds. In fact, as remarked 
earlier, this relation implies the recurrence relation itself. 



4. g-Difference Equations. Our first result is a (/-analogue of Theorem 2.1. 

Theorem 4.1 If w{z) is analytic in the ring q < \z\ < 1 and is continuous on its boundary then 



(4.1) (D q <p n )(z) 



u(C) ~ u(qz) 
C-qz 



+i<t>n(z) 



Cl=l 



u(C) - u(qz) 

C-qz 



Proof. Expand D q <f> n {z) in a series of the n 's. We get 



n-l 



{i- q ){D q 4> n ){z) = I J2m^M0[M0-Mq0M0% 

Break the above integral as a difference of two integrals involving 4> n (C) an d <j>n{qC)i then in the 
second integral replace £ by (,/q. Under such transformation 4>k(C) 1S transformed to 4>k(qC)i since 



|C| = 1. Furthermore (1.4) gives 

(4.2) w((/q) = [l + Cu(C)(l-l/q)MC). 

Therefore 



(l-q)(D q( f> n )(z) 



„ n—l 

4l=ifr L J < 



^n— 1 



y n -i{z) 



., K n— 1 



At : 



4l=l fro < 
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The result now follows from ( |2,9| ). 

We next substitute for (j)* n {z) in (EO) from (pT3|) , if 4> n (0) ^ 0, and establish 



(4.3) {D q( t) n ){z) = AnW^xiz) - B n {z)(j) n {z), 



with 

(4.4) A n {z) 

(4.5) B n (z) 



K n 1 " 



, • K n— 1 

+ *- — —z 



- u(qz) 



(0) J Kl=1 C-qz 



|C|=i C-QZ 



-MO 



w(()d(. 



(o) 



The are the q- analogues of (|2.13| - |2.1q ). Here again we set 



(4.6) L nA = D q + B n {z), 

(4.7) L n>2 = -D q -B n _ 1 (z) + 



The ladder operations are 



A n — (0) 



ZKn-2 



+ 



K n _2^„(0) 



(4.8) L n>1 J> n (z) = A n {z)^ x {z\ L n ^ n ^(z) = ^l^zl^zM^i - 



(/) n (0)K n - 2 Z 



This results in the ^-difference equation 



(4-9) L n , ( Mz) ~ An - liz) ^~ l{ ° )Kn - 



A n {z)' 



Z 4> n {0)K n - 



~4>n{z). 



There is also a ^-analogue of the functional equation ( |3.2| ) which can be found most simply 
by exploiting the third Remark to Theorem (3.1). 



Theorem 4.2 If u(z) is analytic in the annular region q < \z\ < 1 then the following functional 
equation for the coefficients A n (z), B n (z) holds 



(4.10) B n + B 



K n— 1 A n — i 



n-1 



K n 071-1 (0) 



K n -2 Z K n -2 4>n(0) 



A n ~ 



n — 1 u(qz) 1 — q 
qz q q 



n-1 

E 

j=0 L 



B 



K j Aj 

Kj—l Z 
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Proof. Two alternative forms of the second order (/-difference equation are possible, namely 
and the following, 

(4.11) Ln+1,1 ( A Z , s L n+1,2 ) 4>n(z) = — A n+1 (z)(f) n (z). 

These two equations, written out in full are, respectively 

(4.12) E$M*) 

( , . AJqz) , . D a AJz) 



A n (z) ' A n (z) 

A n (qz)A n ^ 1 (z) 

K n (fin— I 

(0) A n (qz)A n ^ 1 (z) 
(0) A n (z) 

A n (qz) A n ^(z)B 

n\ z ) 



D q (j) n (z) 



+ \D q B n (z) - ^f\ Dq A n (z) + ^^B n (z)B n ^(z) 



K n _ 2 A n {z) 



K n (j) n _i(ti) A n (qz) K n -i4> n -i(0) A n (qz)A n _i(z) . 

A n -i(z)B n (z) H — — r } 4> n (z) = 0, 



(0) A n (z) xw " w K n _ 2 n (O) 

and 

(4.13) Z>^„(*) 

( A n (qz) D q A n {z) 



g^ n (z) "" ' ^n(^) 

K n A n (qz) K n+1 <f) n (0) 1 . , 

: A n (gz) + — }> D q (p n (z) 



K n -1 qz Kn-! 4> n +i(0) qz 

D q B n (z) - ^D q A n (z) + #giW^n(.) K " 



(0) A ra (gz)i3 n+ i(z) K n 0n(O) A n (qz)>l n+ i(z) 

K n— 1 0n+l(O) Q K„_l n+ l(O) <?Z 

B n (z) K n+1 4> n (0) A n (qz)\ 

<?Z K n _i0 n+ i(O) J 

A comparison of the coefficients of the first g-difference terms leads to the difference equation 

(4.14) l -B n+1 {z) - B n ^{z) - Kn MZ) + Kn ~ l K ~ 1{Z) 
q K n -\ qz K n _ 2 Z 

- <hM An (z) + -gg- ^W ^cz) = -1. 

Using the results for the first coefficients 

(4.15) B 1 {z) = - u ( qz )-^-M 1 (z), 

(4.16) M^l = - zMl (z), 

K-l 
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with 



M.,7, <"> ,«)* 
JC]=i C - 92 < 



this difference equation can be summed to yield the result in ( 4. 10|) . 



In the next example we will follow the notation and terminology in || and O. The g-shifted 
factorials are 



(a;q) :=l, (a; q) n : = J|(l - aq k x ), n = 1, • • • or oo, 



fc=i 

while the multi-shifted factorials are 



m 



(ai,a 2 ,--- ,a m ;q) n ■= JJ^fc^ 
fc=i 

Define the inner product 



(4.18) (/,<?)=/ /(C)s(C)MC) , ■ 

With respect to this inner product the adjoint of L n i is 



(4.19) (L* n J)(z) = z 2 [q- (l-q)zu(z)]D q f(z) + zf(z) + + u(z)\ f(z), 

provided that w(z) is analytic in q < \z\ < 1 and is continuous on \z\ = 1 and \z\ = q. The proof 
follows from the definition of D q and the fact g(Q = g(l/C), when |£| = 1. Observe that as q — * 1, 
the right-hand side of ( f4.19| ) tends to the right-hand side of Q2,25| ), as expected. 

Example. Consider the Rogers-Szego polynomials {7i n (z\q)} , where 

(4.20) w(z) = — — , H n (z\q) = ^ 



In this case 



n/2 n/2 

(4.21) = -J^—HnWq), 0„(O) 



\/(q;q)n ' y/(q;q)n y/{q\q)n 

It is easy to see that 

(4.22) U (z) = ^- + f^ 
l—o l— a 
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Thus [it(C) — u (q z )]/(C ~ Q z ) is ~~ — l)( z }- A simple calculation gives 

(l_ g n)3/2 Kn _ 1 n _ 1 (z) /■ / a\~7*"7 7T 



(DM(z) = 4>n-l{z) + J / MOKiQCMOir, 

— Q <Pn(U){L-q) J\(\=l l( , 

which simplifies to 

(4.23) (DM (z) = , 



since K n 4>^(qC) — <p n (P)^n{<lO 1S a polynomial of degree n — 1. The functional equation ( 4.23Q can 
be verified independently by direct computation. 



5. Discriminants. Schur [28], [33, §6.71] gave an interesting proof of the Stieltjes-Hilbert 
evaluation of the discriminants of the classical orthogonal polynomials of Hermite, Laguerre, and 
Jacobi. His proof relies on a very clever observation. Let {p n (x)} be a sequence of polynomials 
satisfying a three term recurrence relation 

(5.1) p n (x) = (a n X + b n )p n -i(x) - C n p n - 2 (x), U > 1, 

and the initial conditions 

(5.2) pq(x) = 1; pi(x) = a\x + bt, 

together with the conditions a n _ic n / 0, n > 1. Schur [|2g] observed that 

n n 

(5.3) Upn-xix^) = (-1)»<-1)/* J] <- 2k+1 C k k ~\ 
k=l k=l 

where {xj )TL : 1 < j < n} is the set of zeros of p n {x). 

Let Zj >n be the zeros of 4> n {z). Following Schur, we let 

n 

(5.4) A n = J]> n _i(z jin ). 

5=1 

Lemma 5.1 XYie expression A n is given by 



(5-5) A n = ^n4 n>2, A 1 = l. 
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Proof. It is clear that 



n n—1 



n-l 



(5.6) A n = Yl Y]_(Zj,n ~ Z k ,n-l) = T^TT Y\ <t>n(zk,n-l) 

3=1 k=l ha k=l 



Zk,n—l4 > n—2 (Zk,n-l)- 



The recurrence relation ( |2.5D we find 

K n -2<t>n(ty 



(5.7) (f>n(Zk,n-l) 



«Vi-10ra-l(O) 



Substituting from (|5.7|) into (|5.6[) and applying </> ra _i(0) = K n ~i llj=i ( — we establish the 
two term recurrence relation 



<- 2 [<Mo)] 



n-l 



A n _i, n > 1. 



By direct computation we find Ai = 1, so the above two term recursion implies (|5.5|). 
Examples. For the circular Jacobi polynomials the discriminant is given by 



(5.8) 



—) 

n+a J 



n-l 



(n-l)!(2o+l)„_i 



(«+l)n-l 

while those for the Szego polynomials are 



n/2 g (a+1)? 



^ j!(2a+l) J 



A 



/ a+b 



2rt 



2n-l r 



(5.9) 



\2n+a+b J 
1 



(n-l)!(a+6+l)n-i(a+l/2) w (ft+l/2)n 
(«+^+l)L-l 



i2n-l 



(a+l/2) n (6+l/2), 

2n-2 



Jir^i 1 i!(a+6+l)i(o+l/2)K6+l/2)i 
n-l)!(a+6+l) n _ 1 (a+l/2) n _ 1 (6+l/2) n _ 1 



a — b 



(5.10) 



2n-l+a+6 
x(n-l)!(a+6+l) n _i 



a+^+l)2n-2 



n-l/2 



n?:i 2 («+fe+i) 



Iir=i n(a+6+l),(a+l/2),(6+l/2), 



The resultant of two polynomials f n and g m is 

n 

(5.11) R{fn,g m }=i m l\g m (z J ), 

3=1 



where /„ is as in (|L7|) . Observe that || §100] 
(5.12) D(/„) = (_l)»(»-i)/2 7 -i /;>■ 
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In general let T be a degree reducing operator T, that is (Tf)(x) is a polynomial of exact degree 
Ti — l when / has precise degree n and the leading terms in / and Tf have the same sign. Define 
the generalized discriminant D(f n ,T) by 



(5.i3) D(f n , t) ■.= (-lr^/v^i/n,^} = (-lr^vv-'n™^)' 

3=1 



for f n as in (|L7|) . 



Theorem 5.2 Lei {</> n (z)} 6e orthonormal on the unit circle and assume that T is a linear 
operator such that 

(5.14) T(f> n {z) = A n (z)<j) n -x(z) - B n (z)<j) n {z). 



Let {zk, n ■ 1 < k < n} be the zeros of (j) n (z). Then the generalized discriminant (5. IS) is given by 



I i \n(n— 1)/2 \± (c\X\n— 1 n_1 71 

(5.15) Dfa, T) = Ul JMjjj J] K 2 "Q j4n(zfcn)) „ > o. 

" 1 J=l fc=l 



Proof. Apply ( g3| ), Q and Q. 



In the case of the orthonormal Rogers-Szego polynomials Theorem 4.2 and, (4.21) and (4.23) 
imply the discriminant formula 



n-1 



(5.16) D(<f> n ,D q )=D(<f> n ,q)= [ - q > (q;q) n H 



3=1 



1 



For the Rogers-Szego polynomials we get 



(5.17) D(H n ,q) = (-q) 



-n(n-l)/2 



(1-5) 



n-l 

n 



If one is interested in the limiting case q — ► 1 then we need to rewrite ( |5.17l ) as 



(5.18) D(H n ,q) = (- 



n(n-l)/2/ 1 _ Nn(n-l)/2 



(l-?) n 



n 



(1- 



(?; ?)j 



which shows that D(H n , q) — ► 0, for n > 1, when g — > 1, as expected since 7i n (z|g) — ► (1 + z) n as 
g-1. 
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